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Abstract 

We couple a double-channel imaging technique, allowing for the simultaneous acquisition of 
high-quality and high-resolution intensity and peak emission wavelength profiles [T. Wang and 
G.L. Lippi, Rev. Sci. Instr. 86, 063111 (2015)], to the polarization-resolved analysis of the optical 
emission of a multimode VCSEL. Detailed information on the local wavelength shifts between the 
two polarized components and on the wavelength gradients can be easily gathered. A polarization- 
and position-resolved energy balance can be constructed for each wavelength component, allowing 
in a simple way for a direct analysis of the collected light. Applications to samples, other than 
VCSELs, are suggested. 
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Current technology relies heavily on the high manufacturing quality of light-emitting 
devices, thus on beam quality (spatial and emission wavelength homogeneity) and on polar¬ 
ization control. Among the various light sources, Vertical-Cavity Surface-Emitting Lasers 
(VCSELs) stand out because of their single longitudinal mode operation, small beam diver¬ 
gence, (nearly) circular spot prohle, very low threshold current and ease of integrability. The 
excellent symmetry properties currently obtained in VCSEL cavity structure, contrary to 
the strong asymmetry between TE and TM modes intrinsic to edge emitters, relinquish all 
control on the polarization state of the emitted electromagnetic (e.m.) held - manufacturing 
details and imperfections become therefore of paramount importance and generate complex 
polarization dynamic^®^ where the emission on one linearly polarized component can switch 
to the orthogonal directionP®. The need for polarization characterization and controP^ still 
drives research and development in VCSEL construction. 

In multi-transverse mode VCSELs small frequency splittings, originating from intrinsic 
birefringence, have been identihed as a source of polarization switchin^^ and the minimum 
difference in refractive index between the two components, necessary to stabilize operation in 
one polarization, has been estimated^ to be of the order of 5 x 10“^. Experimental measure¬ 
ments have conhrmed this hndin^^®® and have paved the way for designing polarization- 
stable device^^^^. In the past, numerous techniquehave been devised, including 
interferometric (Fabry-Perot) measurement^, rehectance difference spectroscop}^^^ and 
near-held imagin^^®2ll^ most of them are tailored to a specihc kind of sample or to par¬ 
ticular experimental conditions. Thus, there is still a need for a more general and hexible 
measurement technique. 
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FIG. 1: Experimental setup. The collimator (Thorlabs, C230TM-B) is adjusted to obtain the 
near-held image in the scanning plane of the prohlei® and is preferred here to a microscope 
objective given the minimal distance required to insert the various optical elements (wave plate 
and PBS). The A/2 plate is adjusted to select either polarization component (and the ensemble 
A/2-plate-PBS is removed for unpolarized images). The other key elements - photodetector (high- 
gain UDT-455), the spectrum analyzer (Agilent 86142B) and the digitizing oscilloscope (Agilent 
DSO-X 3024A) - are those already usecP^. Measurements are taken with the spatial resolution 
chosen for the unpolarized scan^2S!_ 

In this letter, we take up this challenge and present a scheme allowing for the simultane¬ 
ous, high-resolution, polarization-resolved cartographic reconstruction of the intensity and 
frequency distribution of emitted e.m. field intensity. The direct imaging of the emission 
frequency, associated to the intensity pattern and to the polarization information, allows for 
an immediate pictorial description of the emission and helps in identifying and quantifying 
the device’s multi-mode features. 

As a test device, we use a multimode, commercial semiconductor laser emitting at A = 
(980 ± 3)nm (Thorlabs VCSEL-980) with maximum output power Pmaxii ~ 10mA) = 
1.85mhF. The measured threshold current is ith ~ 1.3mA and the estimated aperture 
diameter is d ~ 6/im (deduced from the manufacturer’s datasheet). The laser is mounted 
on a temperature-stabilized holder (Thorlabs TCLDM9) with stability better than O.liT 
and is supplied by a commercial source (Thorlabs LDC200VCSEL) with resolution IpA and 
accuracy ±20/iA. The polarization analysis is performed by inserting a A/2-plate, followed 
by a Polarizing Beam Splitter (PBS), in the beam path of a scannning devic^^S! which 
samples the image of the near-held produced by the collimator (Fig. [^. The total average 
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laser power output as a function of injection current and the typical optical spectrum at 
i = bmA are displayed in Fig. which shows: (a) the lasing threshold at ith ~ l.SmA, and 
(b) the emission wavelength centered between 978 and 979nm. 



FIG. 2: Unpolarized input-output laser characteristics (a) and optical spectrum at i = 5mA (b). 

The selected, inexpensive multimode VCSEL is well suited for a detailed cartographic 
demonstration due to its quite complex modal superposition leading to a rich spectral dis¬ 
tribution in the total and polarization-resolved output. Fig. [^a) shows the typical near-held 
intensity distribution emitted by the VCSEL at i = 5mA. A three-fold symmetry is ap¬ 
parent, with a strong localised emission peak in one of the lobes (cf. in particular the 3D 
inset). While the overall symmetry is likely to originate from the superposition of angu¬ 
larly symmetric (Ince-Gauss) mode^, the “hot spot” is probably due to current crowding. 
Fig. I^b) shows the associated, color-coded, spatially resolved peak emission wavelengtfpSl. 
The spectral information complements the one contained in the intensity prohle. Compari¬ 
son between the two panels (Fig. plotted on the same geometrical scale) shows that the 
two upper intensity lobes (panel (a)) share an emission wavelength (A — A^e/ ~ —0.3nm), 
while the bottom one has its peak shifted further towards the blue (A — A^e/ ~ —0.7nm). 
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FIG. 3: Total mean output intensity and wavelength distribution profiles at pumping current 
i = 5mA: (a) intensity distribution profile; (b) spatial wavelength distribution. All spectral 
distributions are plotted - here and in the following - relative to a reference wavelength Xref = 
979.2nm. The color scale in the spectral pictures corresponds to the actual wavelength shift (blue 
representing the shortest wavelength components present in each spectrum). The digit before the 
decimal point in each wavelength scale corresponds to 10“®m (i.e. nm). 

The polarization-resolved input-output laser characteristics are displayed in Fig. and 
show that, as expected, the weaker (Pj_) component has a higher threshold. However, at 
the operating current value, i = 5mA, both components are active, with a power ratio 
P^: Pu^l: 2.5. 


5 








1.4 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

0123456789 



Current (mA) 


FIG. 4: Input-output lasing performance curves for the two Pii and P±_ polarization states. 


The intensity patterns (Fig. [^a,c)) show a polarization selective complementarity in the 
localization of respective maxima and minima, even though they do not lead to a complete 
smoothing of the total output (Fig. [^a)) due to the strongly unequal amount of power 
carried by the two polarizations (Fig. |^. The spectral information, which thanks to the 
spectrometer’s higher sensitivity extends beyond the range displayed by the intensity pro- 
hle (as in Fig. [^b)), shows a much more homogeneous emission in the orthogonal, weak 
component (P^, Fig. |^d)); its emission takes place at one main wavelengths (~ —0.3nm 
- i.e., below Xref) and possesses a basic “rectangular” symmetry. The spectral information 
is much more complex in the principal polarization component (Py, Fig. j^b)) where we 
see the appearence of stronger wavelength shifts in the beam wings with a more elaborate, 
rotationally symmetric structure. 
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FIG. 5: Intensity and wavelength distribution for Py (top) and P± (bottom) polarization states, 
separately, (a) and (c) are the intensity distribution profiles, (b) and (d) display the corresponding 
wavelength distributions. 

The most interesting polarization-related information comes from the spatially resolved 
wavelength difference associated with the two polarized spectra (Fig. [^. The background 
is now set at center scale (O.OOnm) while strong wavelength differences appear over the 
beam’s cross-section (AA||__l = Inm). Strong, isolated, small-sized regions where the emis¬ 
sion wavelength strongly differs from its surrounding appear in the regions bordering the 
beam’s “edges” on the outer border (AA||_j_ > 0) as well as on the inner border (AA||^j_ < 0). 
Their origin may be linked to material inhomogeneities, which would hardly influence the 
laser emission since this points are located in areas where the local intensity is very small 
(compare the coordinates of these points in Fig. and compare to Fig. [^); however, use¬ 
ful information about hner details (device structure, modal superposition, fabrication, etc.) 
may be carried by these spots. Our interest here is simply in detecting and highlighting these 
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regions, irrespective of their physical origin, to demonstrate the sensitivity and resolution of 
our technique. 
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FIG. 6: Spatial distribution of wavelength difference (AA|| _i_ = Ap|| — Apj^) in the overlap region 
for i^i and Pj_ polarization states. 

The cartographic reconstruction well highlights the strong wavelength gradients which 
appear in the regions where the two polarized components switch their local intensity dom¬ 
inance. Interestingly, areas of constant frequency difference between the two polarizations 
appear (e.g., the “keyhole” shape in the center of the beam with wavelength difference 
~ —0.15), while the individual polarization spectra show a more complex pattern 
(Fig. [^,d) and intensity gradients exist in the total (Fig. |^) and polarization-resolved 
(Fig. |^,c) intensity distribution. 

Four “outer” spots (located between 0.6 and 0.7 in the horizontal direction, and 0.4 
and 0.6 in the vertical one in Fig. stand out particularly well and show the presence of 
wavelength gradients up to —~ 0.5nm/fim - estimated from a wavelength difference 
AA||^j_ ~ O.lnm over a spatial interval As 0.2/im (recostructed from the image on the 
basis an estimated diameter d = 6p.m for the VCSEL). Equating the relative wavelength 







variation to the relative refractive index change, we obtain an estimate (for GaAs, n = 3.6) 
for An 0.001. 

Notice that the wavelength resolution depends on the spectrometer. Our current setupP^ 
has shown a peak wavelength sensitivity as low as O.Olnm, i.e., frequency resolution of 3 
GHz (or 13 ^eV in energy units). Better optical spectrum analyzers can further improve 
the performance. This lower limit well matches the requirements for frequency-resolved 
measurements in VGSELs, whose birefringence-induced splitting has been found to rangd^ 
between 3GHz and 22GHz. It is therefore envisageable to use this technique to obtain 
spatially-resolved spectral information and reconstruct possible localised structural defects 
inducing birefringence in VGSELs with this technique. 

The cartographic demonstration has been based on a multimode VGSEL, but the tech¬ 
nique is suitable for many different samples. For instance, photoluminescence measurements, 
routinely conducted in semiconductor^ and thin hlmd^^® , but also for monitoring the 
quality of large, centimeter-sized boules of semiconductor materiald^S, could beneht from 
the high spatial resolution, coupled to polarization and spectral information offered here. 

Before concluding, from Figs. [^b,d) we estimate the contribution of the different fre¬ 
quency components to the whole beam: Fig. [^distinctly shows the presence of different sets 
of discrete emission frequencies for the two polarized components. This type of information, 
which provides statistical information on the relative fractions of the beam emitting at dif¬ 
ferent wavelengths, can be straightforwardly completed by an energy balance (not shown): 
adding up the intensity in each pixel, as a function of peak emission wavelength, one can end 
up with information equivalent to a peak spectral reconstruction, with the added advantage 
of the knowledge of the position from which each contribution is derived. In other words, 
from this plot one can quantify the amount of energy carried at each frequency, in a fashion 
similar to an optical spectrum (e.g.. Fig. [^b)) where in addition one can also correlate the 
information to the geometrical location of each contribution. 
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FIG. 7: Histogram of the number of “pixels” (i.e. measured points in the sampled matrix) for each 
peak wavelength for Pm (black) and P± (red) polarization components. 


In conclusion, we have demonstrated a cartographic technique which provides 
polarization- and spatially-resolved intensity and wavelength distributions of the emitted 
light with a resolution which can reach 3GHz. A demonstration has been given specihcally 
using a multimode VCSEL, but the technique can be applied to virtually any light-emitting 
source. 
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